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Abstract. Nambu-Goto model is investigated by using the canonical Monte Carlo simulation technique on 
dynamically triangulated surfaces of spherical topology. We find that the model has four distinct phases; 
crumpled, branched-polymer, linear, and tubular. The linear phase and the tubular phase appear to be 
separated by a first-order transition. It is also found that there is no long-range two-dimensional order in 
the model. In fact, no smooth surface can be seen in the whole region of the curvature modulus a, which 
is the coefficient of the deficit angle term in the Hamiltonian. The bending energy, which is not included 
in the Hamiltonian, remains large even at sufficiently large a in the tubular phase. On the other hand, the 
surface is spontaneously compactified into a one-dimensional smooth curve in the linear phase; one of the 
two degrees of freedom shrinks, and the other degree of freedom remains along the curve. Moreover, we 
find that the rotational symmetry of the model is spontaneously broken in the tubular phase just as in the 
same model on the fixed connectivity surfaces. 



PACS. 64.60.-i General studies of phase transitions 
tronic - 87.16.Dg Membranes, bilayers, and vesicles 



1 Introduction 



Triangulated surfaces are one of the basic models to in- 
vestigate the physics of biological membranes and that 
of strings [Tl[51[31|311Sl[nil7] . Surface models can exhibit a 
variety of shapes because of the two-dimensional nature. 
A well-known model is the one of Helfrich, Polyakov and 
Kleinert (HPK) [8,9, 10 . Studies have focused on the phase 
structure of the surface model of HPK 11, 12, 13 , 14, 15, 16, 
fTZlfTSl , and it was reported that the model has the smooth 
phase and the crumpled (or folded) phase on the fixed- 
connectivity/fluid surfaces [l9l[2f)l[2T1[2^[28l[24l[25l[26l[27l[28l 
l29..30.,3L 32.3a..34.,35j. 

A branched-polymer phase and some linear structures, 
which are not always smooth, can be seen in a fluid model 
on dynamically triangulated surfaces [36j . A surface model 
of Nambu-Goto also has a variety of phases; the crumpled 
phase, the tubular phase, and the smooth phase, even on 
the fixed connectivity surfaces of spherical topology when 
the Hamiltonian includes a deficit angle term [37], which 
is an intrinsic curvature energy |i38lli39[|40[|41| . 

However, little attention has been given to the linear 
structure (such as a smooth curve) in the surface models. 
Therefore, we study in this paper the Nambu-Goto model 
|42j with the deficit angle term on dynamically triangu- 
lated surfaces of spherical topology. 



68.60.-p Physical properties of thin films, nonelec- 



The Nambu-Goto surface model is well known as an ill- 
defined one whenever the Hamiltonian includes not only 
no additional term but also the standard bending energy 
term [43]. However, we have already confirmed that the 
Nambu-Goto surface model changes to a well-defined one 
on the fixed connectivity surfaces when the deficit angle 
term is included in the Hamiltonian |37j as well as a cer- 
tain bending energy term |34| . Then, it is also interesting 
to see whether or not such well-definedness remains un- 
affected on the dynamically triangulated fluid surfaces in 
the whole range of the curvature coefflcient a(> 0) of the 
deficit angle term. From the standard argument for fluid 
membranes, we expect that the thermal fluctuations be- 
come large on the fluid surfaces. For this reason, we think 
it is worthwhile to study whether the model is well-defined 
or not on the fiuid surfaces. 

The purpose of this study is to understand the phase 
structure of the surface model of Nambu-Goto with the 
deficit angle term on dynamically triangulated surfaces 
of spherical topology. It is also aimed to see whether the 
model is well defined on the fluid surfaces. 

It will be shown in this paper that the model is well- 
deflned and has four distinct phases; the crumpled phase, 
the branched-polymer phase, the linear phase, and the 
tubular phase, which consecutively appear as the coeffi- 
cient a increases from a sufficiently small value (a = 10) 
to a sufficiently large one (a~ 1 x 10"*). We will find that 
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the crumpled phase, the branched-polymer phase, and the 
linear phase are connected by higher-order transitions at 
relatively small a, and that the linear phase and the tubu- 
lar phase are connected by a first-order transition at rela- 
tively large a. 

A remarkable result is the appearance of the linear 
phase, which has not yet been found in the surface mod- 
els of HPK and in the fixed connectivity Nambu-Goto sur- 
face model. In the linear phase of the model in this paper, 
one of the two dimensions of the surface shrinks around a 
curve, and the other dimension remains along the curve. 
Moreover, a spontaneous breakdown of the rotational sym- 
metry is also found in the tubular phase, where the surface 
spans and extends along a one-dimensional straight line. 
We have no smooth-phase (or two-dimensional order) in 
the model even at sufficiently large a in contrast to the 
same model on the fixed connectivity surfaces [37] . 

2 Model 

The partition function Z of the discrete Nambu-Goto sur- 
face model is defined on a triangulated surface and is given 

by 

N 

T •' i=l 

S{X,T)^Si+aS3, (1) 

A i 

where / flfc^i '^-^i ^^le 3iV-dimensional integrations, and 
^-j- denotes the sum over all possible triangulations T. 
The Hamiltonian S is given by the linear combination of 
the area energy term 5*1 and the deficit angle term 6*3 
such that S{X,T) = Si+aSa. S{X,T) denotes that S 
depends on the embedding X and the triangulation T of 
the triangulated surface of spherical topology. In the area 
energy Si in Eq.((T]), A/^. is the area of the triangle A. The 
symbol Si in S3 is the sum of the angles of vertices meeting 
at the vertex i. 

The center of the surface is fixed in the integration 
in Eq.ll]) to remove the translational zero mode. Z{a) of 
Eq.([T]) denotes that the model is dependent on the curva- 
ture coefficient a. It should also be noted that the Hamil- 
tonian S{X, T) is defined only with intrinsic variables of 
the surface. 

The term 5*1 in Eq. ([1]) is a discretization of the origi- 
nal Nambu-Goto action defined by 5* = / dPx^, where 
g is the determinant of the first fundamental form on 
the surface X swept out by strings. The surface X is lo- 
cally understood as a mapping from a two-dimensional 
parameter space into R'^. The mapping X is not always 
injective, because the surface is allowed to self-intersect 
in string models. Therefore, the model can be called the 
phantom surface model as long as it is considered as a 
model for membranes in R'^. Note also that the Polyakov 
string model can be obtained from the original Nambu- 
Goto model [44] . 



The deficit angle term 5*3 is closely related to the in- 
tegration measure dXi We have to remind our- 
selves of that dXi can be replaced by the weighted mea- 
sure dXiqf, where qi is the co-ordination number of the 
vertex i, and a is considered to be a = 3/2. Considering 
that Qi is a volume weight of the vertex i, we assume a as 
an arbitrary number. Moreover, the co-ordination number 
Qi can be replaced by the vertex angle 6i such that 

n,dX,q^ -> n,dX, exp(a J2 log ^0 ■ (2) 

i 

The constant term ^ ^ log 2t: is included to normalize S3 in 
Eq.(IT]) so that = when 6i — 2tt at the vertices. Thus, 
we have the expression S'3 in Eq.([T]). Note that S'3 = 
is satisfied not only on the flat surface but also on the 
cylindrical (or tubular) surfaces. 

The unit of physical quantities is as follows: The length 
unit a in the model can arbitrarily be fixed because of 
the scale invariant property of the partition function. As 
a consequence, the string tension coefficient X in S = 
A5i-l-aS'3 can be fixed to A=l, because A has the unit of 
kT/a'^, where k is the Boltzmann constant and T is the 
temperature. The coefficient a has the unit of kT. 



3 Monte Carlo technique 

The icosahedron is used to construct the triangulated lat- 
tices for the starting configuration of the Monte Carlo 
(MC) simulations. The edges of the icosahedron are di- 
vided into £ pieces of uniform length, then we have a tri- 
angulated lattice of size iV = 10^^ -I- 2. 

The canonical Metropolis technique is used to update 
X, which is the three-dimensional vertex position of the 
triangulated surface. The position Xi is shifted such that 
X'^—Xi + AXi, where AXi is randomly chosen in a small 
sphere. The new position X^ is accepted with the prob- 
ability Min[l,exp(-Z\5')], where Z\5' = S'(new) - S'(old). 
The radius of the small sphere for AXi is chosen at the 
beginning of the simulations to maintain 35 '--^ 65 % ac- 
ceptance rate; almost all MC simulations are performed 
on about 50 % acceptance rate. 

The variable T is summed over by the bond flip tech- 
nique, which destroys the uniform lattice structure of the 
starting configurations constructed from the icosahedron. 
The vertices freely diffuse over the surface and, there- 
fore, such surfaces are called the fluid random surface. 
A sequential number labeling the vertices and another se- 
quential number labeling the bonds become at random, 
because the bond flips change the pairing between a bond 
and the corresponding two vertices. By using the sequen- 
tial number of bonds, a bond is randomly selected to 
be flipped, and the flip is accepted with the probability 
Min[l,exp {—AS)]; this is one update of T. N updates of 
X and N updates of T are consecutively performed and 
make one MCS (Monte Carlo sweep). The acceptance rate 
rr of the bond flip is not fixed a priori. We have rr = 0.65 
at a = 10 in the crumpled phase, and it increases with 
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(a) a=10(Af = 2892) 



(b) a = 50(iV = 2892) 



increasing a and is almost independent of a at a > 30, 
where 7-7- = 0.53 ~ 0.55. 

The lower bound 10~^ is imposed on the area of tri- 
angles in the updates of X and T, where Aq is the mean 
area of the triangles computed at every 500 MCS, and Aq 
remains constant due to the relation Si/N = 1.5. We re- 
mark that the areas are almost free from the lower bound, 
because the areas of almost all triangles are larger than 
lO^^Ao throughout the MC simulations. No constraint is 
imposed on the bond length. A sequence of random num- 
bers called Mersene- Twister [57] is used to update X and 
T in the simulations. 



4 Results 

(a) Q = 10(iV=1442) 



(b) a = 50(iV = 1442) 




(c) a = 2000(iV = 1442) (d) a = 3000(iV = 1442) 

Fig. 1. Snapshots of the A'^ — 1442 surfaces obtained at (a) 
a = 10 (crumpled phase), (b) a = 50 (branched-polymer phase), 
(c) a = 2000 (linear phase), and (d) a = 3000 (tubular phase). 
Figures (a),(b), and (d) were drawn in the same scale, which 
is different from that of (c). The thickness of the surface in (c) 
is drawn many times larger than the original one, which is too 
thin to draw. 



Snapshots of the iV==1442 surfaces are shown in Figs. 
[TJa), [Ijb), [Ijc), and [ijd), where the surfaces were re- 
spectively obtained at a = 10 (crumpled phase), a = 50 
(branched-polymer phase), a — 2000 (linear phase), and 
a = 3000 (tubular phase). Figures [DJa), [l^b), and [T^d) 
were drawn in the same scale, which is different from that 
of Fig[TJc). The thickness of the linear surface in Figlljc) 
is drawn many times larger than the original one, which 
is too thin to draw. The reason why we call the surface 
in Figmb) as the branched-polymer surface is because of 
its shape. Although the surface in Figdfb) appears to be 
almost linear, we can see that some parts of the surface 
are branched. 




(c) a = 5000(Af = 2892) 

Fig. 2. Snapshots of the iV 



(d) a = 6000(iV = 2892) 

= 2892 surfaces obtained at (a) 



q: = 10 (crumpled phase), (b) a = 50 (branched-polymer phase), 
(c) a = 5000 (linear phase), and (d) q = 6QQ0 (tubular phase). 
Figures (a),(b), and (d) were drawn in the same scale, which 
is different from that of (c). The thickness of the surface in (c) 
is drawn many times larger than the original one, which is too 
thin to draw. 



Figures EJa), EJb), EJc), and M,<i) show snapshots of 
the N = 2892 surface, which were obtained respectively 
at a = 10 (crumpled phase), a = 50 (branched-polymer 
phase), a — 5000 (linear phase), and a — 6000 (tubular 
phase). Figures m a), [H^b), and O^d) were drawn in the 
same scale, while the scale of Figl2ljc) is different from 
that of the other three figures. The thickness of the linear 
surface in Fig|2{c) is drawn many times larger than the 
original one, which is too thin to draw. 

The mean square size X^ is defined by 



X' 



(3) 



where X is the center of the surface. 

Figure [3fa) shows X^ obtained at relatively small a 
close to the phase boundary between the crumpled phase 
and the branched-polymer phase, where the surface size 
is iV = 812, N = 1442 and N = 2892. In FigElJb), we 
plot X^ obtained at relatively large a close to the phase 
boundary between the linear phase and the tubular phase, 
where three different sizes; N = 812, N = 1442 and N = 
2882, were also assumed. The curves in FiglJl^a) show that 
the surface size continuously increases with increasing a. 
On the contrary, X^ shown in FiglSjb) abruptly changes 
and then, it clearly indicates a discontinuous transition 
between the linear phase and the tubular phase. 

The dashed lines drawn vertically in both of the figures 
represent the phase boundaries of the N = 1442 surface. 
The position of the phase boundary in Figl3{b) is imme- 
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Fig. 3. (a) obtained at relatively small a close to the 
boundary between the crumpled phase and the branched- 
polymer phase, and (b) obtained at relatively large a close 
to the boundary between the linear phase and the tubular 
phase. The dashed lines drawn vertically in both of the fig- 
ures represent the phase boundaries of the N = 1442 surface. 

diately found between the linear phase and the tubular 
phase, because X"^ discontinuously changes at the bound- 
ary. On the contrary, the phase boundary between the 
crumpled phase and the branched-polymer phase in Fig[3fa) 
is unclear. Therefore, the position of it was roughly deter- 
mined by viewing the snapshots of surfaces at the bound- 
ary region of a. 

The branched-polymer phase and the linear phase are 
also ambiguously separated; they are expected to connect 
smoothly to each other, because continuously changes 
at the corresponding region of a as can be seen in Figs|3l[a) 
and[3jb). The position of the phase boundary of the iV = 
1442 surface is at a = 80 ~ 100, which was also expected 
by viewing the snapshots of surfaces. 

The convergence speed of MC simulations is very low 
in the linear phase close to the tubular phase. About 6 x 
10^ 8 X 10^ MCS were done for the thernialization at 
a = 3000 5000 on the TV = 2892 surface. One reason 
of such low convergence speed in the linear phase seems 
due to the large phase-space volume in R'^. The crumpled 
surfaces in the crumpled phase occupy a relatively small 
region in R'^, while long, thin, and string- like surfaces in 
the linear phase can extend to a large space in R'^. Then, 
the canonical and local update procedure of X is very 
time-consuming for such string-like surfaces to have the 
equilibrium configurations as long as the surface obeys 
the Hamiltonian of short ranged interactions. 

We comment on the Hausdorff dimension H , which is 
defined by the relation 

^2 oc N^/" (4) 

at sufficiently large A^. Figures HJa) and [Hb) are log-log 
plots of X^ against N . The data X"^ in Fig|4ja) were ob- 
tained at a = 10 in the crumpled phase and at a = 50 in 
the branched-polymer phase, and the data X^ in FigdJ^b) 
were obtained at a where X^ is at its maximum in the 
linear phase and at a in the tubular phase close to the 
linear phase. Thus, we have 

H^w = 4.32 ± 0.26, = 3.30 ± 0.44, 




Fig. 4. (a) Log-log plots of against A'^ obtained at a = 10 
and a = 50, and (b) log-log plots of X^ against N obtained at 
a where X^ is at its meiximum in the linear phase and at a in 
the tubular phase close to the linear phase. 



iJlin = 1.52 ± 0.02, iJtub = 1.56 ± 0.06. (5) 

It is quite natural that the result Ha=w — 4.32(26) is 
larger than the physical bound _ff = 3 in the crumpled 
phase, because the surfaces are allowed to self-intersect 
and completely crumpled. The result Ha^o ~ 3.30(44) 
is larger than the value H = 2, which is specific to the 
branched-polymer surfaces. This indicates that the branched- 
polymer phase of the model is not exactly the branched- 
polymer phase. However, we call the surface obtained at 
that region of a as the branched-polymer surface because 
of its surface-shape such as shown in Figs[Hb) and[2l^b), 
as stated above. 

We find that iJii„, which was obtained in the linear 
phase close to the tubular phase, is almost identical to 
-fftub, which was obtained in the tubular phase close to 
the linear phase. We understand that both results iJun = 
1.52(2) and i?tub = 1-56(6) are larger than the value H—1 
of a one-dimensional straight line of constant density of 
vertices and less than the value H = 2 oi the branched 
polymer surfaces, and moreover, Hn^ and i?tub are slightly 
larger than 7J= 1.22(33) in the tubular phase of the fixed 
connectivity surface model [37| . 

The size of the surface can also be measured by the 
maximum linear size D, which can be approximated as 
follows: Firstly, find the vertex / that has the maximum 
distance from X the center of surface, and secondly, find 
the vertex J that has the maxi mum distanc e D from the 
vertex I. Then, we have D = yJ{Xi - Xj^. Note that D 
is identical to the diameter when the surface is a sphere. 

Figuresl^a) and[5l[b) show D of the surfaces at a close 
to the phase boundaries. The size D changes discontinu- 
ously at the transition point in Figin^b). We can also un- 
derstand from Figini^b) that the string-like surfaces curve 
and entangle themselves to occupy a spherical region in 
R3 in the linear phase, since D(2892)/i:»(1442)~ 1.37(lin) 
is greater than 1.26~ (2892/1442)1/3 and less than 2.01~ 
2892/1442. The value 1.26 for the ratio L>(2892)/i:>(1442) 
is given by assuming that the string-like surface, whose 
length is proportional to N , forms a spherical region in R'^. 
If the region has a constant density of vertices, then the di- 
ameter of the region is proportional to N^/^ . In this case. 
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Fig. 5. The maximum hnear size D obtained at (a) relatively 
small a close to the boundary between the crumpled phase 
and the branched-polymer phase and at (b) relatively large ct 
close to the boundary between the linear phase and the tubular 
phase. The dashed lines drawn vertically in both of the figures 
represent the phase boundaries of the A'^ = 1442 surface. 



change can be seen in 5*2 /Nb ■ On the contrary, 5*2 /Nb in 
FigElJb) appears to change discontinuously at the transi- 
tion point between the linear phase and the tubular phase. 
This indicates that the transition is of first order, although 
the value of 5*2 /Nb is relatively large even in the tubular 
phase. The large value of 82/ Nb implies that the surface 
is not smooth even at a = 8000 on the iV=1442 surface. 

It is interesting to see whether the smooth surface ap- 
pears at sufficiently large a. Therefore, we performed MC 
simulations at a = 1 x 10*, a = 2 x 10*, a = 3 x 10*, 
a = 5 X 10*, a = 1 X 10^ on the iV = 1442 surface. As a 
consequence, we found that the obtained 82/ Nb are al- 
most identical to those in the tubular phase in Figl^b). 
This implies that the smooth surface can be seen only 
a,t a —^ 00. Thus, we conclude that no smooth phase can 
be seen in the model on fluid surfaces in contrast to the 
same model on the fixed connectivity surfaces, where the 
smooth phase can be seen at sufficiently large but finite 
a. 



the corresponding Hausdorff dimension is H = 3. On the 
other hand, the value 2.01 for the ratio i:)(2892)/i:)(1442) 
is given by assuming that the string-like surface spans a 
one-dimensional straight line of constant density. The cor- 
responding Hausdorff dimension is iJ = 1 in this case. 

We also find from FigOb) that D(2892)/i:»(1442) ~ 
1.53(tub) in the tubular phase close to the Hnear phase. 
The value 1.53(tub) is greater than 1.42 ~ ^2892/1442 
and less than 2.01 ~ 2892/1442. The value 1.42 (2.01) for 
the ratio i5(2892)/£)(1442) is given by assuming that the 
surface forms a tubular surface, whose density of vertices is 
proportional to N^^^ (N^) per unit-length, and therefore, 
the corresponding Hausdorff dimension is H = 2 {H = l). 
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Fig. 6. The bending energy S2/NB obtained at (a) relatively 
small Q and at (b) relatively large a, where A''^ is the total 
number of bonds. The bending energy 52 is not included in 
the Hamiltonian. 



The bending energy S2/NB is plotted in FigsElJa) and 
inUb), where Nb(^ 3N — 6) is the total number of bonds, 
and S2 is given by S'2 = X]y(l-nrnj), where Y,ij denotes 
the summation over unit normal vectors and rij of the 
triangles i and j sharing a common bond. Although S2 is 
not included in the Hamiltonian, it can reflect how smooth 
the surface is. We find from FiglHI^a) that no discontinuous 
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Fig. 7. The area energy Si/N obtained at (a) relatively small 
a and at (b) relatively large a. 



In order to see that the model is well-defined, we plot 
Si/N versus a in FigslJlJa) andjT^b). We find from the 
figures that Si/N satisfies the expected relation Si/N — 
3{N — l)/2iV~3/2. If the relation was violated, then the 
equilibrium statistical condition should not be expected 
in the model. In fact, we are unable to see the relation 
Si/N = 3/2 in an ill-defined model such as the one defined 
only by the term Si in Eq.([T]). Thus, the results shown 
FigsIZUa) and [7fb) indicate that the model in this paper 
is well-defined. 

The value of the deficit angle term per vertex S3/N is 
plotted in FigslH^a) andlUJb). If the transition is of first or- 
der, one can see a discontinuity in S3/N. However, S3/N 
appears to change smoothly against a at the phase bound- 
aries. We hardly see a discontinuous change in S3/N even 
at the discontinuous transition point in FiglSfb). One rea- 
son of this seems come from the fact that the coefficient a 
is very large while S3 is very small at the transition point. 
In fact, the Hamiltonian is able to have a finite jump even 
when 5*3 has a small jump because of the large value of a. 
If we have many simulation data at the phase boundary 
between the linear phase and the tubular phase, a discon- 
tinuity can be seen in S3/N. However, we must note that 
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00 
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crumpled 



(a) 

Fig. 8. The deficit angle term S3/N obtained at (a) relatively 
small a and at (b) relatively large a. 



the possibility of continuous transition is not completely 
eliminated. 



Fig. 10. (a) The histogram h{L) of the bond length L of tri- 
angles at the phase boundary between the crumpled phase and 
the branched-polymer phase, and (a) those at the phase bound- 
ary between the linear phase and the tubular phase. The sam- 
ples of the bond length L were obtained at every 1000 MCS in 
the final 2 x lO'^ MCS. 
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Fig. 9. The maximum bond length Lmax obtained at (a) rela- 
tively small a, and at (b) relatively large a. Lmax were obtained 
in the final 2 x lO'^ MCS at each a. 



The maximum bond length -Lmax is expected to change 
depending on a. In contrast to the standard Gaussian 
bond potential, the area energy Si gives a constraint only 
on the area of triangles and no constraint on the bond 
length in the Nambu-Goto model. Consequently, long and 
thin triangles appear in the equilibrium configurations and 
form such linear surfaces. Therefore, we have to check that 
L-max < D] the maximum bond length is less than the lin- 
ear size D except in the crumpled phase, where D can be 
comparable to Lmax because of the crumpled nature. In 
order to see this, we plot L^ax versus a in FigslSIa) and 
EJb). These i^ax were obtained in the final 2 x 10^ MCS 
at each a. We find from Figsl9]and[5]that L^ax is compa- 
rable to D in the crumpled phase as expected, and that 
■^^max is quite larger than D in the linear phase. Moreover, 
imax is still larger than D in the tubular phase. In fact, 
D = 79, Lniax = 55 at a = 3000 on the N = 1442 surface, 
and D = 121, L,nax = 66 at a = 6000 on the N = 2892 
surface. 

Not only the surface shape but also the surface struc- 
ture can be influenced by the phase transitions and changes 
depending on the value of a. Figures fWa) andfTUfb) show 
the distribution (or histogram) h{L) of the bond length 



L obtained at a close to the phase boundaries. The his- 
togram h{L) was obtained during the final 2 x 10^ MCS 
at each a, and h{L) is normalized such that the graph 
h{L) and the horizontal axis enclose a constant area. We 
find from FigfTUTa) that h{L) at a = 10 (the crumpled 
phase) is different from h{L) at a = 50 (the branched- 
polymer phase). This fact indicates that the lattice struc- 
ture in the crumpled phase is different from that of the 
branched-polymer phase, although h[L) is expected to 
change smoothly from one phase to the other phase. 

It is also found from Fig fTOT b) that the histogram h{L) 
in the linear phase (a = 3000) is different from the one 
in the tubular phase [a = 6000). We can expect that 
the shape of h{L) discontinuously changes because of the 
discontinuous nature of the transition, in fact, it is ob- 
vious that h{L) at one phase boundary in Fig fTUf a) is 
quite different from h{L) at another phase boundary in 
Fig flpy b). However, h[L) shown in Fig fTOT b) appears to 
vary smoothly from one phase to the other. The reason 
of this seems because the surfaces are composed of ob- 
long triangles both in the linear phase and in the tubular 
phase. We must recall that h{L) in the fixed connectivity 
model discontinuously changes at the transition point be- 
tween the tubular phase and the smooth phase [37]; the 
reason of the discontinuous change of h{L) is because the 
triangles in the tubular phase are oblong while those in 
the smooth phase are almost regular. 

We remark that the histogram h{A) of the area A of 
triangles is not influenced by the phase transitions in con- 
trast to h{L), just as in the flxed connectivity model in 
[37] . In fact, we checked that no difference can be seen 
in the histograms h(A) obtained in the crumpled phase, 
in the branched-polymer phase, in the linear phase, and 
in the tubular phase on the surfaces of size N = 812, 
TV = 1442, and iV = 2892. 



5 Summary and conclusion 

A surface model deflned by the Nambu-Goto Hamilto- 
nian is investigated by MC simulations on dynamically 
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triangulated surfaces of size up to iV = 2892. The Hamil- 
tonian includes a deficit angle term. The Nambu-Goto 
surface model is well known as an ill-defined one if the 
Hamiltonian includes not only no additional term but also 
the standard bending energy term. It was also reported 
that the deficit angle term makes the Nambu-Goto sur- 
face model well-defined on the fixed connectivity surfaces 
of spherical topology, and the phase structure of the model 
was clarified [37] . 

In this paper, we aimed at showing whether the model 
is well defined or not on the fluid surfaces and how the 
phase structure changes when the surface changes from 
the fixed connectivity one to the fluid one. The fluid sur- 
face model in this paper has a single parameter a just like 
in the fixed connectivity model in [37 . The only difference 
between the model in [37j and the one in this paper is that 
the model in [37] is defined on the fixed connectivity sur- 
face while the model in this paper is defined on the fluid 
surface. 

We found that the fluid model in this paper is well de- 
fined and has four distinct phases; the crumpled phase, the 
branched-polymer phase, the linear phase, and the tubu- 
lar phase. The first three are smoothly connected, and the 
last two are connected by a first-order phase transition. 
We should note that the possibility of higher-order transi- 
tion between the linear phase and the tubular phase was 
not completely eliminated. There is no smooth phase in 
the whole range of «(> 0). 

One remarkable result is that the model has the linear 
phase, where the surface shrinks to such a one-dimensional 
curve. One of the two-dimensions of the surface seems to 
be spontaneously compactified, and the remaining dimen- 
sion survives along the curve. A spontaneous breakdown 
of the rotational symmetry is also found in the tubular 
phase, where the surface spans a tubular surface along a 
one-dimensional straight line. 

Moreover, the phase transition is reflected not only in 
the surface shape but also in the surface structure. In fact, 
the histogram h(V) of the bond length varies against a 
and, h,(L) in one phase is different from those in the other 
phases, although seems smoothly varies even at the 
transition point between the linear phase and the tubular 
phase. This point is in sharp contrast to that of the fixed 
connectivity model, where h(V) changes discontinuously 
at the transition point between the tubular phase and the 
smooth phase [37]. On the contrary, the histogram h{A) of 
the area A is not influenced by phase transitions as in the 
fixed connectivity model, and no difference can be seen in 
h(^A) in four different phases. 

This work is supported in part by a Grant-in- Aid for 
Scientific Research from Japan Society for the Promotion 
of Science. 
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